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Abstract 

n 

We applied the approach of K.-Y. Yang, T. M. Rice and F.-Ch. Zhang [1] (YRZ) to analyze the 
high resolution angular resolved photo-emission spectroscopy (ARPES) data in BiSCO obtained 
recently at Brookhaven. In the YRZ ansatz a constant RVB gap is assumed which leads to Luttinger 
zeros along the AFBZ and four Fermi pockets centered on the nodal directions. We relax the 
assumption of a constant RVB gap function, treating it as a Ising order parameter accompanied by 
thermal fluctuations. If these thermal fluctuations are very strong leading to strictly short range 
correlations in the spatial dependence of the RVB gap, then the reconstruction of the Fermi surface 
into pockets will not survive. We examined the intermediate case of critical fluctuations leading 
to a power law falloff of the RVB gap correlations. To this end we followed the analysis recently 
developed by two of us to treat the effect of the power law correlations in the antiferromagnetic 
2-dimensional xy model on the single electron Greens function. The partial truncation of the Fermi 
surface to form pockets survives in the presence of power law correlations. The linewidth of the 
quasiparticle peaks increases with the exponent of the power law correlations. If this exponent is set 
at the value of critical fluctuations in the 2-dimensional Ising model a relatively small linewidth is 
obtained. If this exponent is doubled, a linewidth comparable to the values found in recent ARPES 
experiments on underdoped BiSCO (Tc = 65 K) at T = 140K was obtained. The anisotropic 
suppression of the quasiparticle peak around the Fermi pockets hides the back side of the pockets 
closest to the AFBZ to give essentially Fermi arcs seen experimentally. 

PACS numbers: 71.10.Pm, 72.80.Sk 



I. INTRODUCTION 

The underdoped regime of the cuprates remains the greatest mystery of condensed matter 
physics. Its resolution may have very profound influence on our understanding of strongly 
correlated phenomena. For several years after their discovery it was believed that the 
cuprates had a large hole like Fermi surface (FS) which area is proportional to 1 — x with 
X being concentration of the holes. This was celebrated as a triumph of the Luttinger the- 
orem which apparently held even in such strongly correlated systems. Gradually more and 
more conflicting experimental evidence have emerged. Now it is overwhelmingly obvious 
that in their normal state the underdoped cuprates have a FS consisting of either finite arcs 
or pockets centered around nodes of the Brillouin zone. The evidence for arcs comes from 
ARPE_S exp.i,ne„t. fll. and *e evidence fo. pockets co^es f.o™ bo.h .ece,. ARPES 
data k| and most conclusively from quantum oscillation measurements in high magnetic 
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fields |5|,|6|. 

As far as the theory is concerned, there are two camps - traditional and radical ones. 
The traditional camp attempts to explain the phenomenon within the familiar paradigm 
of physics of weakly interacting systems. The radical camp argues that the doped Mott 
insulators such as cuprates display essentially nonperturbative features which do not occur 
in weakly interacting systems. The first question to ask is whether this phenomenon is a 
manifestation of some new physics characteristic for strongly correlated systems or can be 
understood within the traditional paradigm. 

In the traditional scenario the pockets originate from the Fermi surface reconstruction 
occurring in the presence of some density wave order parameter. This order parameter is 
likely to compete with the superconductivity and several candidates have been suggested. 
However, with the exception of LBCO at special doping no static ordering leading to increase 
of the unit cell has been observed. This difficulty is explained away with a statement that 
for FS reconstruction to occur it is sufficient to have a fluctuating order with a sufficiently 
large correlation length. Detailed calculations in support of this scenario based on the so- 
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called spin-fermion model can be found in the paper by Chubukov and Morr 

more recent paper by Sedrakyan and Chubukov |8|. As the fluctuating order parameter the 

authors have taken commensurate spin density wave (SDW). This choice leads to certain 

difficulties. 



Perhaps, the most serious one comes from the fact that fluctuating SDW assumes 
predominance of soft modes in the spin excitation spectrum. However, the pseudogap 
phenomenon is associated primarily with the decrease of magnetic susceptibility which 
speaks against such a scenario. What is important is that such a decrease has been 
observed in clean stoichiometric cuprates such as YBa2Cu408 [9] and HgBa2Cu04+5 
(lO|,[ll|. This excludes an interpretation of the spin gap as disorder phenomenon. 
Static CDW or SDW order would also lead to a splitting of NMR lines. Meanwhile, 
such effects have not been seen in the aforementioned stoichiometric cuprates. To the 
contrary, the cited experiments show narrow NMR lines on different oxygen sites and 
a clear evidence of the spin gap in the temperature dependence of the Knight shift and 
the Gaussian spin-echo decay rate. It is true that due to the technical problems related 
to the surface preparation ARPES and STS experiments on thee materials have not 
been done. This leaves a possibility to argue that there are no Fermi pockets there. 
This is, however, a narrow escape route since in all other respects these materials are 
not different from other cuprates. 

The other difficulty is related to the fact that the choice of commensurate SDW places 
the centers of the pockets squarely at (±7r/2, ±7r/2) points of the Brillouin zone (Fig. 
1). This seems to be incompatible with the recent high resolution ARPES measure- 
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The latter difficulty does not go away even if one invokes incommensurate spin density 
wave suggested in [24], [12]: 

S = 3f?e[NieiQ^'- + N2e^^^'"], 

Qi = (7r-^», Q2 = (7r,7r-5) (1) 

Although such SDW is naturally conjectured from the neutron experiments, which 



observe incommensurate low energy peaks in dynamical magnetic susceptibility 13| , in 
such choice would not do the job. Due to the Umklapp processes there are not four, but 
eight hot spots on the bare (large) FS connected by the wave vectors (±7r±(5, ±it±6). 
This arrangement again produces the electron spectral weight symmetric with respect 
to the magnetic Brillouin zone boundary. 




FIG. 1: Formation of a Fermi pocket in the mean field Fermi surface reconstruction caused by 
the SDW ordering. The bare Fermi surface (thin solid line) is not nested. The dashed line is 
the magnetic Brillouin zone boundary. Two subbands with i = 1,2 are connected by the anti- 
ferromagnetic wave vector Q = (vr,7r). 



Nontraditional scenarios directly appeal to strongly correlated physics and suggest that 
in strongly correlated systems such as doped Mott insulators FS may be small incorporating 
only holes and this does not require any ordering. Dzyaloshinskii reminded us that the 
Luttinger theorem, in fact, does not relate the particle density to the FS volume 1^ (see 
also [15]). The latter relation exists only in weakly correlated systems where the electron self 
energy S(ci;,k) does not have singularities at zero frequency. In strongly correlated systems 
where S(a; = 0, k) goes to infinity at certain surface in momentum space, the volume inside 
of this surface contributes to the particle density. A Mott insulator is an obvious example of 
an insulating state with a half filled Brillouin zone which seemingly violates the Luttinger 
theorem (in its traditional, or rather, according to Dzyaloshinskii, confused understanding). 

This idea was illustrated by Essler and Tsvelik [16] who considered a strongly correlated 
model which allowed a controlled approximation. This was a model of Hubbard chains 
coupled by a long range tunneling (the inverse tunneling radius k served as a small parameter 
of the theory). Strong intrachain correlations were taken into account non-perturbatively and 
the interchain tunneling and exchange interaction processes were considered in a controlled 
fashion by RPA and perturbation theory in k. The salient feature of this model was existence 
of a small FS in the form of hole and electrons pockets. These pockets existed in a finite 
temperature range; below a certain temperature the model undergoes a magnetic ordering. 



'ield (CDMF) approach; states with small FS have been found (see, for instance, 



The main purpose of this model was to illustrate the principle. Namely the lack of the 
relation between the carrier density and the volume of the FS. To apply these ideas to the 
cuprates one needs to study more realistic models. This was done using Cluster Dynamical 

VIean 

l3,[l8|,ll9|). 

In principle, the radical scenario does not assume an extreme proximity to ordered states, 
though the model of Essler and Tsvelik has gapless collective excitations. Konik, Rice and 
Tsvelik (KRT) [20] obtained pockets in a model of coupled ladders which does not have 
gapless modes though it is fair to say that this model has two sites per unit cell and hence 
the translational invariance was already broken. Based on these results Yang, Rice and 
Zhang (YRZ) [ij] conjectured a phenomenological expression for the single particle Green's 
function which turned out to be highly successful in describing the ARPES and other data 



in underdoped cuprates, e. g. see |2l|.|22|: 



GKk)=[(. + iO)-.(k)-^-J|i^]-, (2) 

where A(k) = (cos k^ — cos ky)A and ^o = io(cos k^ + cos ky) and e(k) = ti(cos k^ + cos ky) + 
t2 cos k^ cos ky — fi describes the bare band. The function A(k) is the main input in the 
theory and requires explanation. In the KRT model the parent state is the so-called d-Mott 
state of two-leg ladder which has no local order parameter. Instead, at T = it has a 
topological order and A is the quasiparticle gap. It is not apriori clear what is the status of 
this function in two-dimensional case. One can imagine that it becomes a fluctuating order 
parameter and this is the option we explore in this paper. 

An expression similar to ([2]) has recently been obtained by Qi and Sachdev J23i] on the 



basis on the theory of algebraic charge liquid (ACL) suggested by Moon and Sachdev [24 1. 
The presence of an ACL places this theory into the radical camp. The difficulty is that the 
ACL state has gapless collective nonmagnetic excitations which presence would presumably 
be acutely felt tn the thetmodyna^ic. As a mattet of fact, the theoty of the pocket 
formation suggested in [23i] is not that different from [161]. In both cases the parent high 
energy state is fractionalized; the excitations are spinless particles with charge e and neutral 
particles with spin S=l/2. Their bound states are quasiparticles with a small FS. An 
essential difference comes from the status of the parent state. For Hubbard chain this state 
is well established through the exact solution of this model. It is less certain with ACL. The 



only example of such a state has been found in the model which physical parameters are 
dramatically different from those for the cuprates [25J. Ironically, this is a situation typical 
for the cuprate physics: those models which allow rigorous treatment do not apply and those 
models which may apply do not allow rigorous treatment. The other difference comes from 
the influence of gapless excitations. In [l6| the interaction of quasiparticles with gapless 
collective modes strongly renormalizes the self energy resulting in a marginal Fermi liquid 



state. In [23| the interaction between gapless gauge field excitations and quasiparticles is 



weak and the state is essentially a Fermi liquid. 

II. EXPERIMENTAL DATA AND FITTING PROCEDURE 

In this paper we report an attempt to fit the ARPES data for BiSCCO using a modified 
form of YRZ spectral density. The ARPES spectra shown in this paper were recorded on 
beamline U13UB at the NSLS using a Scienta SES2002 electron spectrometer. The samples 
are cleaved and measured in a UHV chamber, where the base pressure is maintained at 
5 X 10^^^ Torr. Spectra are recorded in the form of spectral intensity plots which represent 
a map of the intensity as a function of binding energy and angle of emission or momentum. 
The energy and angular resolution used in these studies was 15 meV and 0.1° respectively. 
The experimental curves presented by dashed lines on Figs. |3], JH El have been obtained 
from the data by the deconvolution procedure described in [^]. This procedure removes 
a substantial part of the broadening caused by the experimental resolution effects. The 
spectral intensity may be sampled at a given binding energy as a function of momentum, also 
known as momentum distribution curves (MDC). The Bi-2212 crystals were grown by the 
floating zone method. To get the under-doped BSCCO single crystals, the optimally doped 
material was annealed under vacuum at 550 °C for 48 hours. The T^ of underdoped single 
crystals is 65 K, which was determined using SQUID magnetometry. The measurements 
were performed at fixed temperature T =140K. 

The logic of our approach is as follows. Our goal is to check whether fluctuation effects 
are compatible with the formation of FS pockets. To this end we perform a kind of stress 
test on the YRZ and similar radical scenarios. In accordance with these scenarios we will 



follow the approach taken in 26|] and consider the spectral function of quasiparticles moving 



in a slowly varying order parameter field. Further we assume that the fluctuations are as 
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FIG. 2: (a) The hole Pocket can be approximated as an ehipse with seniiaxes a and h. The rotated 
system of coordinates (A:^, /cm), with A;_l measured along the nodal direction is used to present the 
data, (b) The zoom of the tip of the pocket region which contains the cuts 2 and 3. Points k*^^^ 
and k*^^-* on the pocket boundary are characterized by maximal ARPES intensity in MDC data 
along cut 2 and 3. Sufficiently close to the tip of the pocket, k*-^-* it can be approximated by the 
hyperbolic curve asymptotically approaching two lines symmetric with respect to the kn direction 
and crossing at the point k^*^'. The hyperbola touches the pocket and represent accurately the 
front side of the pocket. The back side of the pocket gets shifted away from the line k^ = 0. This 
shift of the ghost FS is of little consequence as the ARPES intensity is very small at the back 
side of the pocket. The angle (/> formed by the bare Fermi surface with the direction along fcii is 
determined by fitting the experimental data taken along cut 2 and cut 3. Vectors vi and V2 are 
the velocities at the real and ghost FS respectively. 



weak as possible in the absence of long range order. This means fluctuations with longest 
possible correlations length, i.e. critical ones. At finite temperatures critical fluctuations 
are also classical. So we assume that the system in question has a wide critical region where 
the order parameter fluctuations have power law correlations. In this way we consciously 
underestimate the effects of fluctuations. However, if we manage to fit the data, this would 
mean that these radical approaches have a chance. If not a more radical scenario is needed. 
As far as the order parameter is concerned, we believe that the considerations described 
above exclude most of the explanations with a possible exception of [23\. However, the 
latter theory also points to weak fluctuation effects. 



To fit the experimental data we adopt the approach of Ref. 26j to treat classical fluctu- 
ations in YRZ theory. In Ref. 26|] the model of electrons in two subbands of opposite spins 



with dispersion relations ei(k) = e(k), and 62 (k) = e(k + Q) coupled by the antiferromag- 
netic order parameter has been studied, see Fig. [H The effective Hamiltonian in the two 
subband space describing this situation can be written as 

H = ei(k)^^ + e2(k)i^ + A(k)^ + A*(k)^ , (3) 

where cxj are Pauli matrices acting in the two subbands space, and cr± = cr^ dz iay. This 
Hamiltonian with constant order parameter A(k) corresponds to the mean field solution. In 
this limit of no fluctuations the Fermi surface is reconstructed and the Green function reads 

G„,(..k)=[(. + iO)-.,(k)-^;J|g^]-\ (4) 



To study the effects of fluctuations, it has been assumed in [26|] that the spin magnetization 
either i) has a strong easy plane anisotropy and lies in the XF-plane or ii) has a strong easy 
axis anisotropy and is directed along the z axis. Most importantly the regime of critical 
fluctuation has been studied in both cases. In the first case the fluctuations are critical in 
the whole range of temperatures below a BKT transition. While in the second case they are 
critical only at the transition. The self-energy calculated to second order in perturbation 
theory in A yields for the Green function 

where d is the scaling dimension of A, and £ is the ultraviolet scale of the order of the lattice 
spacing. The new energy scale entering the Green function in Eq. ([5]) is T^ = A(A£)°'/^^~'^). 
The self energy in Eq. (J5]) is regular at the mass shell, u = ei(k). However at higher orders 
it contains singular terms at both the mass shell, u = ei(k) and at the ghost mass shell 
u = 62 (k). These singularities were systematically resummed in Ref. [26i]. We focus the 
discussion on the result obtained close to the mass shell, where most of the spectral weight 
accumulates, and comment on the behavior at the ghost Fermi surface later. In the former 
case the result reads 

i_ I 2dt'^T^{2~2d)A^ . / G^^ \ 2di£'^'^T^{2 - 2d)T{l - 2d)A^ 
^ + {-i{u - e2(k)))4-4d 2 ^(^^_,2(k)J+ (-i(a; - e2(k)))5-6'^ ' 

(6) 

where G2 is given by Eq. ([5]). Equation (|6]) is vahd in the limit \uj — ei(k)| < |w — e2(k)|, 
and Tk <^\uj — e2(k)|. 



As we have said, the expression ([6]) is quite general. If the bare Fermi surface is considered 
as the first subband , ei(k) = e(k), and we take the magnetic zone boundary for the second 
one, 62 (k) = ^o(k) Eq. 1^ is obtained in the mean field, Eq. (^. With this correspondence, 
Eqs. (E]), dS]) become 

i_ 1 2de^'^r\2-2d)A^ 1 / G^^ \ 2dii^'^r\2 - 2rf)r(l - 2d)A^ 

(8) 
These equations have a meaning of the phenomenological YRZ function modified by fluctu- 
ations. In this modification we treat A in ([2]) as an Ising-like order parameter and assume 
that its fluctuations are classical and critical. The choice of Ising order parameter dictates 
the value of the scaling dimension of A as d = 1/8. However in the absence of full mi- 
croscopic derivation of Eq. (j2]) the scaling dimension d should be really considered as a 
fitting parameter with fluctuations enhanced for larger vales of d. In fact, we have found 
that d = 1/4 more accurately describes our data. We additionally assume that fluctuations 
remain critical in a wide temperature range as is the case for the XY model below the BKT 
transition. This assumption probably introduces a bias towards order thus making the spec- 
ial lines narrower than they must be. We remind the reader that the scenario described in 



23|, where quasiparticles have ~ w^ + T^ decay rate, is likely to lead to even sharper peaks. 
Below we introduce simplifying approximations to the dispersion to facilitate the com- 
parison of Eqs. ([7]), ([8]) with the experimental data. Close to the tip of the pocket the bare 
dispersion e(k) can be linearized, e(k) ^ vi(k — k^*^^). Here the point k^*^) is the crossing of 
the bare Fermi surface e(k) = with the line passing through the two tips of the pocket, 
see Fig. [2l In what follows it is convenient to choose k.^^^ as the origin of coordinate frame 
in the Brillouin zone with axes fcy and k± directed parallel and perpendicular to the mag- 
netic zone boundary respectively, see Fig. [2j In the above coordinate frame the velocity 
Vi = {v sin 0, —V cos 0) is specified by the angle it forms with the diagonal direction, see 
Fig. O For simplicity we also modify the line of zeros of the Green function in Eq. ([21), 
^o(k) — 7- V2k with V2 = {vsm(f),vcos(f)). We stress that although this modification leads 
to a deformation of the back side of the pocket which is discussed later, it does not change 
appreciably the behavior of the Green function at the mass shell which is the focus of current 
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discussion. In the linear approximation the mean field, pocket is approximated by the hy- 
perbolic curve asymptotically approaching lines perpendicular to velocities vi^2, see Fig. |2l 
The requirement that the hyperbola touches the pocket and has the same curvature at the 
tip imposes a relation between the distance |k*^^^ — k^^^l and the angle (p. The latter is the 
only free parameter for fixed d which has to be determined by the fitting procedure described 
in App. [Al 

We express the result ([8]) in the universal form G = T^^G(k, a)) through dimensionless 
momentum, k = vh/Tx and energy u = u/Tk, with 



' -\{ijj — k±_ cos — k\\ sin 0))*^""^'^ \uj — k± cos — k\\ sin < 



G-i =G^i + -— = ^ ^, TT^^G^^ In 



2diT^{2 - 2d)T{l - 2d) 
(— i(ci} — k± cos — k\\ sin 0))5~6d ' 



where 



G-i = c^ + ^x cos - fc|| sin + , ., , '^^\ ^'^^ . ,,,,„,, . (10) 

(— i^w — /e_|_ COS0 — /e|| sm0)j^ ^"^ 

To compare our results with experiment we plot the spectral density, A{u,k\\,k±) = 
—l/7TlmG{u,k\\,k±) for fixed k\\ as a function of k± corresponding to MDCs taken along 
the cut across the pocket. We first consider the FS, w = 0. In Fig. |3] we show the fit 
to the experimental data shown by dashed lines with the theory for the scaling dimension 
d = 0.125. This value would correspond to a transition region of two-dimensional Ising 
model, or to the BKT transition in the XY model. The theoretical expression ([8]) used in 
the comparison describes the spectral density at the real FS. In Fig. [3] it is plotted in the 
interval k± > which includes the real FS, or more generally the mass shell, u = e(k). 
Both the data and the fit are presented for cut 2 and cut 3 taken parallel to the nodal 
(F-X) direction, see Fig. [2] with two different values of k^\. Positions of both peaks are fitted 
simultaneously. The angle and the ratio k±/k_i_ = T^/v used to plot Figs. [3], HI [5l [6] are 
obtained through the approximate procedure described in App. |Al This gives an energy 
scale Tk = 54meV. We notice that the obtained energy scale is of the order of pseudogap, 
Tk ~ T* as expected. Although the position of the peaks is captured by the fit their width, 
is smaller than in the actual data. The main point is that one can clearly see from Fig. [3] that 
the FS survives critical fiuctuations. For the negative energy u = — 45meV corresponding to 
the hole states well below the FS the data is shown in Fig. H] (dashed line). The theoretical 
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fits are plotted according to the same Eq. ([8]) with previously obtained scale Tk- Here 
again the peak position is reproduced quite closely, but theory yields peaks which are too 
narrow. We argue that other sources of broadening such as Landau damping mentioned 
earlier might be responsible for the extra width of the signal. To see if the FS survives in 
the presence of stronger fiuctuations we consider larger values of the scaling dimension d. In 
the absence of detailed microscopic derivation of Eq. ([2]) this parameter is considered as free. 
In Fig. |5]we show the fit to the experimental data taken at the FS, u = with d = 0.25. 
In this case the agreement is improved as the higher values of d correspond to a stronger 
fiuctuations, and in result give a broader MDC. The same fitting procedure (see App. |X]) 
yields a slightly different value for the energy scale, T^ ~ 67meV. In Fig. E] we show the 
fit with identical parameters but for a negative energy, u = — 45meV. At this energy the 
fit is much less accurate, but it still refiects the basic observation that the FS withstand 
the critical fiuctuations. To summarize, critical fiuctuations slightly modify the shape of 
the pockets, introduce a finite smearing as is refiected in the shape of MDCs. But most 
importantly their effect is not strong enough as to restore the FS topology, and pockets 
survive. 

Until now we have focused on the region k± > which contains the FS. In the present 
paragraph we discuss briefiy the fate of the back side of the pockets lying in the comple- 
mentary interval, k± < 0. The critical fiuctuations tend to suppress the ghost FS, so that 
the spectral weight mostly accumulates at the front side of the pocket or close to the real 
FS, leading to Fermi arcs. We still discuss the Green function close to the ghost Fermi 
surface, u ~ ^o(k) here for the sake of completeness. Equation ([8]) is not valid anymore in 



the interval k± < and has to be replaced. In the Ref. |26| close to the ghost FS, u ~ e2(k), 
we have obtained the following result 



G 



1 A2(k)£2V^r(l-2rf) 



a;-ei(k) (a;-ei(k))^ 



u;-e2(k)-f- 



iA\k)f^T{l - 2d) 
(-i(a;-ei(k)))i-2d 



(11) 



The result in Eq. f lT2|) was obtained in the following way. In the second order of perturbation 
theory in A, Eq. ([5]), the Green function acquires a correction singular at the ghost FS, 
u = 62 (k). Higher orders corrections are even more singular. Equation ( ITTl) is result of 
careful resummation of the most singular terms at a; = e2(k). It is therefore applicable in 
the limit jw — e2(k)| < |u; — ei(k)|, and T^ ^ |u; — ei(k)|. In exact similarity to the previously 
considered case of the FS in order to adopt the result (ITT|) to our present purposes one has 
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to make a substitution ei(k) = e(k), and e2(k) = ^o(k). The resulting expression reads 



1 A\k)f%^^'^T{l-2d) 

u-e{k) ^ (a;-e(k))2 



w-eo(k)4 



iA^{k)f'^T{l-2d) 



(12) 



;-i(a; - e(k)))i-2'^ 

In the mean field limit, d = 0, the expression (TT2|) goes over to the Eq. ([2]) as expected. At 
finite d it expresses the effect of critical fluctuations on the YRZ ansatz at the ghost FS in 
much the same way as Eq. ([8]) describes the modification of Eq. ([2]) by critical fiuctuations 
at the real FS. The applicability conditions of Eq. ( IT2|) break down sufficiently close to the 
tip of the pocket, k ^ k*^^^ see Fig. [2l and our results are expected to be less reliable in this 
region. Indeed, the better theoretical fits are obtained for the cut 3 taken further away from 
the tip of the pocket than for the cut 2 in Figs. [3l HI [5] and [61 

In order to enable the comparison of Eq. ( IT2l) with experiment we repeat the steps leading 
to Eq. (jH]). Namely we rewrite Eq. ( IT2l) in the form G = T^^G(k,u}) with dimensionless 
momentum, k = vk/T^ and energy u = uj/Tk, 

-_ 1 ei^'='r(l-2rf) 

^ —~ f ', f — : — r + 



Lo + k± cos (j) — kii sin (f) ' {lo + k± cos (j) — k\\ sin (f))"^ 

1+2(1 

■ (13) 



X 



'- T ^ T ■ j^ ir(l-2rf) 

uj — k± cos — k\\ sm (p + 



-i(a; + A;_L cos(/) — ^11 sin0))^ '^^ 

The energy scale Tk the angle in Eq. ( !T3|) are identical to those in Eqs. ([8]), (|9]). We 
emphasize that in general the spectral density at the ghost Fermi surface tends to be sup- 
pressed by fiuctuations with most of the spectral weight found in the region A;^ > 0, close 
to the Fermi surface. We use the form of Eq. f lT5]) to obtain the spectra for k^ < and 
the results are included in Figs. El IH and El The spectra for the two regions A;^ > and 
k± < mostly match each other. The discontinuity at fc^ = is smaller for cut 3 than 
for cut 2, see Fig. |2l This is to be expected since as discussed above both Eq. ([8|) used for 
A;^ > and Eq. (fT2|) used for A;^ < are better satisfied away form the tip of the pocket. 

III. CONCLUSIONS 

In this paper we have examined the robustness of radical theories which have been pro- 
posed to describe the anomalous properties of the pseudogap state. We label these as radical 
because they are not derived from a broken symmetry order. In the YRZ case a phenomeno- 
logical form for the single particle Greens function is put forward in analogy to the form 
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FIG. 3: Theoretical fits shown in sohd hne, to experimental data shown in dashed line taken along 
the Cut 2 (red) and Cut 3 (black) in Fig. [2] at zero energy with d = 0.125. The expression Eq. ([9]) 
(Eq. (fTBj) ) is used to plot the theoretical curves for A;_|_ > {k± < 0). The fit in a subinterval of 
negative k± is not shown where the result Eq ()13p is not applicable. 



that was derived for a two dimensional array of 2-leg Hubbard ladders near half-filling. In 
the ACL case 23|] the key ingredient is the choice of the quantum disordering of a antiferro- 
magnetic state in the presence of doped holes. In both cases the Fermi surface is partially 
truncated and a qualitative description of key experiments of the pseudogap state is ob- 
tained. An important open issue is the effect of thermal fluctuations and the robustness of 
the key properties in their presence. In the YRZ ansatz a constant RVB gap is assumed 
which leads to Luttinger zeros along the AFBZ and four Fermi pockets centered on the nodal 
directions. We relax the assumption of a constant RVB gap function, treating it as a Ising 
order parameter accompanied by thermal fluctuations. If these thermal fluctuations are very 
strong leading to strictly short range correlations in the spatial dependence of the RVB gap, 
then the reconstruction of the Fermi surface into pockets will not survive. In this paper we 
examined the intermediate case of critical fluctuations leading to a power law falloff of the 
RVB gap correlations. To this end we followed the analysis recently developed by two of 
us to treat the effect of the power law correlations in the antiferromagnetic 2-dimensional 
xy model on the single electron electron Greens function. The main conclusion is that the 
partial truncation of the Fermi surface to form pockets survives in the presence of power 
law correlations. The linewidth of the quasiparticle peaks increases with the exponent of 
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E = -ibmeV; d = 0.125.Ciit2, 3 




FIG. 4: Theoretical fits shown in sohd hne, to experimental data shown in dashed line taken along 
the Cut 2 (red) and Cut 3 (black) in Fig. [2] at the energy E = -45meV with d = 0.125. The 
expression Eq. Q (Eq. (USD) is used to plot the theoretical curves for A:^ > {k^ < 0). The fit in 
a subinterval of negative k± is not shown where the result Eq (J13p is not applicable. 

the power law correlations. If this exponent is set at the value of critical fluctuations in 
the 2- dimensional Ising model a relatively small linewidth is obtained. If this exponent is 
doubled, a linewidth comparable to the values found in recent ARPES experiments on un- 
derdoped BiSCO (Tc = 65 K) at T = 140K was obtained. The anisotropic suppression of 
the quasiparticle peak around the Fermi pockets hides the back side of the pockets closest 
to the AFBZ to give essentially Fermi arcs seen experimentally. We conclude that the main 
features of the Fermi surface reconstruction in a radical theory such as the YRZ ansatz, are 
robust against power law thermal fluctuations. 

We finally comment on the situation when there is no true order and the opposite case of 
an ordered state. In the absence of true order the fluctuations are characterized by the finite 
correlation length. The pockets don't survive and the FS remains large and not truncated. 
In the case of critical fluctuations studied in the present paper the self energy includes both 
real and imaginary part. The real part modifies the dispersion of the quasi-particles leading 
to a FS reconstruction. The imaginary part produces finite width of the spectral function. 
The imaginary part is smaller than the real part by a factor of d. In result, for d < 1 the 
FS pockets are well defined. The situation is different in the absence of order: while the 
fluctuations still produce a flnite broadening, there is no real part leading to a FS truncation. 
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E = 0; d = 0.25.C)it2,3 




0.05 0.1 0.15 0.2 



FIG. 5: Theoretical fits shown in sohd hne, to experimental data shown in dashed line taken along 
the Cut 2 (red) and Cut 3 (black) in Fig. [2] at zero energy with d = 0.25. The expression Eq. ([9]) 
(Eq. (fTBj) ) is used to plot the theoretical curves for /c_|_ > (A;_l < 0). The fit in a subinterval of 
negative /c_l is not shown where the result Eq (J13p is not applicable. 

For that reason the pockets are destroyed in the absence of order. 

The influence of thermal fluctuations on the spectral density in the ordered state has 
been studied recently in Ref. [8|. The authors have employed the eikonal approximation 
and summed up the most singular corrections in the parameter T / J\ loge|, where J is the 
exchange interaction in the Heisenberg model, and e is the anisotropy parameter of a quasi- 
2D system. In the presence of long-range order the pockets have been found. In addition, 
some signatures of pockets remain even without long-range order. The fluctuations are 
strong enough in this case to redistribute the spectral weight evenly between the real and 
the ghost FS. Our situation is intermediate in the sense that we have only quasi long range 
order, but the fluctuations have infinite correlation length. Our conclusion is that this 
situation is rather similar to the ordered case, i. e. the FS still has small pockets. The 



results of the Ref. 



indicate that we have considered the strongest fluctuations possible 



consistent with small FS topology. 
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E = -45meV: d = 0.25. Cut2,?, 



-is 




0.05 0.1 0.15 0.2 



FIG. 6: Theoretical fits shown in sohd hne, to experimental data shown in dashed line taken along 
the Cut 2 (red) and Cut 3 (black) in Fig. [2] at the energy E = -45meV with d = 0.25. The 
expression Eq. (^ (Eq. (fT^ ) is used to plot the theoretical curves for k± > {k± < 0). The fit in 
a subinterval of negative k± is not shown where the result Eq (|13p is not applicable. 

Acknowledgments 

We are grateful to A. V. Chubukov, R. M. Konik, I. A. Zaliznyak for valuable discus- 
sions. We acknowledge support of the Center for Emergent Superconductivity, an Energy 
Frontier Research Center supported by the US DOE, Office of Basic Energy Sciences. M. 
K. acknowledges support from BNL LDRD under Grant No. 08-002. 

Appendix A: Fitting procedure 



In this appendix we describe the fitting procedure used to estimate the characteristic 
energy scale T^- In our scheme the only fitting parameter is the critical dimension d. For 
a given d we approximate the Fermi surface by the hyperbolic curve touching the pocket 
at its tip, k = 'k^^\ see Fig. [21 obtained for the linearized bare dispersion. The bare Fermi 
surface is characterized by the angle (p it forms with the magnetic Brillouin zone boundary. 
The pocket in turn is approximated by an ellipse with major and minor semiaxes a and b 
respectively. 



In order to fix we find k]_ , and k]_ by solving the equation ReGg ^(w = 0, k]^ \ /cn ) = 0, 
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i = 2,3 numerically for different values of (f). From the Fig. |2]we deduce ku = ku + 0.11a, 
and /cm = /cm + 0.36a. Using the geometrical relation, 1 = k\ [a/lP') tan^ and placing the 
tip of the hyperbolic curve at ki ^ l/sin0 obtained form Re(j2^^(a; = 0, A;^ = 0, A;||) = 

— (n\ — /■ox 

at the mean field c? = 0, the ratio k\_ /k\_ can be found. The latter becomes equal to the 
experimental value k]_ /k]_ at ~ 0.2l7r for d = 0.125. In view of various approximations 
made in our procedure our estimate for requires a fine tuning and the value ~ 0.237r for 
d = 0.125, and ^ 0.227r for d = 0.25 better fit the data. The ratio kf'^ /kf^ = kf Ikf^ = 
Tk/v gives Tk = 54 meV for d = 0.125 and T^ = 67 meV for and d = 0.25 respectively. 
The scale Tk is of the order of magnitude of the pseudogap energy T*. 
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